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ON  PREDICTING  THE  THRESHOLD  CURRENT  DENSITY  OF 
QUANTUM  WELL  DIODE  LASERS 


By 

Gabgyu  Lim 
May,  1992 


Chairman:  Peter  S.  Zory 

Major  Department:  Electrical  Engineering 

The  threshold  current  density  of  GaAs/AlGaAs  graded 
refractive  index  (GRIN),  separate  confinement  heterostructure  (SCH), 
single  quantum  well  (SQW)  laser  diodes  was  measured  at  25°C  and 
compared  with  predicted  Jth  using  published  CZR  (Chinn,  Zory, 
Reisinger)  model.  The  results  indicated  that  the  existing  CZR  model 
for  QW  lasers  was  found  to  be  deficient. 

To  sort  out  the  problem,  we  chose  two  cases,  thick  40nm  well 
(bulk)  and  thin  7.5  nm  QW.  In  thick  wells,  all  the  complicated 
phenomena  related  to  the  QW  vanish  and  only  carrier  scattering  (CS) 
is  expected  to  be  fundamental.  We  investigated  the  effect  of  CS  on  Jth 
using  standard  bulk  theory.  In  thin  7.5  nm  QW  where  quantum  size 
and  nonradiative  effects  are  important,  we  studied  the  effect  of  CS  on 
Jth  using  new  valence  band  mixing  (VBM)  QW  theory. 


vi 


A reasonable  fit  to  the  thick  QW(bulk)  experimental  data  at 
both  25°C  and  125°C  could  be  obtained  using  strict  k-selection  bulk 
theory  and  a knowledge  of  the  carrier  injection  efficiency^  jn).  In 

making  the  fit,  CS  effect  is  not  necessary,  quite  contrary  , to  Yamada's 
results.  This  finding  suggested  the  possibility  that  the  effect  of  CS  on 
Jth  might  not  be  important  even  in  thin  QW  lasers.  This  was  verified 
with  VBM  QW  theory  combined  with  a knowledge  of  the  carrier 
injection  efficiency^ jn).  Combination  of  the  effects  by  valence  band 
mixing  and  measured  injection  efficiency  fit  the  7.5  nm  QW 
experimental  data,  leaving  no  room  for  the  CS  contribution.  Leakage 
current  density  estimated  in  terms  of  measured  injection  efficiency 
(riin)  is  significant,  comparable  to  the  current  density  contribution  by 
CS  in  the  CZR  model.  This  does  not  necessarily  rule  out  the  possibility 
that  experimental  spectra  can,  in  principle,  be  obtained  which  are 
subject  to  the  strict  k-selection  rule. 

It  is  concluded  that  good  predictions  of  threshold  current 
density  at  any  temperature  can  be  made  for  QW  lasers  provided  one 
uses  QW  theory  developed  in  this  work  and  one  knows  how  to 
predict  the  temperature  dependence  of  rijn. 


vii 


CHAPTER  1 
INTRODUCTION 


The  idealized  light  output  power  versus  injection  current 
relationship  for  semiconductor  lasers  is  a straight  line  as  shown  in 
Fig.  1-1.  For  high  performance,  the  threshold  current  (Ith)  should  be 
small  and  slope  efficiency  of  the  line  big.  It  is  well  known  that  the 
quantum  well  (QW)  design,  characterized  by  a very  thin  active  layer, 
is  best  for  high  performance. 

Various  aerospace  applications  have  been  envisoned  in  which 
high  performance  QW  diode  lasers  are  the  key  elements.  However, 
many  of  these  applications  require  that  the  devices  operate  reliably 
at  temperatures  as  high  as  125°C.  QW  lasers  can  perform  better  than 
thick  active  layer  double  heterostructure  lasers  at  high  temperature, 
but  predicting  their  performance  is  much  more  complicated  than 
predicting  that  for  DH  lasers  because  of  their  complex  internal 
structure.  The  goal  of  this  study  was  to  integrate  the  results  of 
experimental  and  theoretical  work  to  determine  a model  capable  of 
predicting  threshold  current  density  (Jth)  at  any  temperature  for  QW 
lasers. 

Many  publications  have  been  devoted  to  theoretical  gain  vs 
threshold  current  calculations  through  the  years.  The  Stern  model 
[Ster73]  was  largely  accepted  as  the  correct  model  for  such 
calculations  in  the  70s,  although  it  predicted  Jth  values  not  precisely 
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Fig.  1-1.  Idealized  light/current  characteristic  of  a semiconductor 
laser  showing  threshold  current  Ith  and  linear  slope. 


3 


the  same  as  those  measured  by  the  best  lasers  of  that  time  [Thom80]. 
In  the  early  ’80s,  Yamada  pointed  out  that  Stern's  spectral  gain 
function  was  too  optimistic  because  it  ignored  carrier  scattering  (CS). 
By  convolving  the  spectral  gain  function  with  a Lorentzian 
broadening  function  (CS  time  of  about  10‘13  sec),  Yamada  and 
Ishiguro  were  able  to  fit  their  threshold  data  for  DH  lasers  [Yama81]. 
Since  then,  most  Jth  modeling  work  for  both  DH  and  QW  lasers  have 
incorporated  the  CS  concept  [Asada84], 

The  above  modelling  works  have  been  based  on  the  parabolic 
model,  which  assumes  a simple  parabolic  band  for  both  conduction 
and  valence  band  structures  and  a constant  optical  matrix  element 
for  the  transitions  anywhere  in  the  bands.  More  recently,  Chinn, 
Zory  and  Reisinger  published  a detailed  parabolic  model(  = CZR 
model)  [Chin88,  Chin89],  where  valence  band  nonparabolicity  was 
included  for  the  light  hole  states.  Other  important  features  included 
in  this  model  were  intraband  scattering,  gain  anisotropy, 
nonradiative  processes  such  as  leakage  of  the  carriers  over  the 
confining  potential  barriers  [Case78,  Dutt81],  and  Auger 
recombination  [Take85]. 

One  success  of  the  CZR  model  was  to  explain,  semi- 
quantitatively,  the  anomalous  temperature  behaviour  of  Jth  for  GaAs, 
thin  QW  diode  lasers  observed  by  Zory  et  al.  [Zory86a].  It  was  also 
used  to  predict  Jth  and  its  temperature  dependence  (specified  by 
characteristic  temperature  T0)  for  30/70  GaAs/AlGaAs  graded 
refractive  index  (GRIN),  separate  confinement  heterostructure  (SCH), 
single  quantum  well  (SQW)  lasers  near  room  temperature  [Chin89]. 
In  order  to  verify  the  validity  of  the  CZR  model  over  a wide  range  of 
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active  layer  thickness  for  30/70  GaAs/AlGaAs  GRIN-  SCH-SQW  laser 
diodes,  we  compared  experimental  Jth  and  T0  values  with  the  above 
predictions.  The  results  showed  that  the  CZR  model  can  predict  Jth 
reasonably  well  for  some  cases.  However,  it  underestimates  Jth  in 
thin  QW  lasers,  and  overestimates  Jth  substantially  for  thick  QW 
lasers.  The  trend  of  T0  values  agrees  pretty  well  with  experimental 
values  but  they  are  overestimated  because  Jth  is  not  predicted 

properly.  The  possible  reasons  why  the  CZR  model  predictions  do  not 
agree  with  experiment  are  as  follows:  (i)  CZR  model  breaks  down  for 
thick  well  (bulk),  (ii)  Lineshape  function  used,  which  was  assumed  to 
be  independent  of  well  thickness,  temperature,  carrier  energy  and 
carrier  density,  is  not  appropriate,  (iii)  Nonradiative  current  might 

not  be  estimated  correctly,  (iv)  CZR  model  does  not  treat  the  valence 
band  correctly.  Fundamental  studies  [Cola87,  Ahn88]  have  shown 
that  the  valence  band  structure  of  quantum  wells  deviates 
significantly  from  that  of  the  parabolic  band  due  to  valence  band 

mixing  (VBM)  as  the  QW  width  decreases.  VBM  gives  nonparabolic 
bands  which  affect  both  the  subband  energies  and  the  density  of 
states  functions.  The  density  of  states  functions  are  found  to  differ 
considerably  from  the  simple  step-like  shape  computed  from  the 

band  edge  effective  masses.  The  band  mixing  also  results  in  large 
differences  in  the  momentum  vector  dependence  of  the  matrix 
elements,  and  moreover,  some  of  the  transitions  which  were 
previously  assumed  to  be  forbidden  become  partially  allowed.  The 
quantum  well  gain  calculated  by  the  k.p  method  and  k-selection 
rules  is  considerably  different  than  those  from  the  conventional 
methods  with  parabolic  bands  and  with  matrix  elements  which  have 
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constant  total  magnitude.  Although  the  effects  of  the  band-mixing  on 
the  band  structure  and  gain  of  a QW  have  been  studied  [Cola87, 
Ahn88,  Zory92],  very  little  work  has  been  done  on  the  comparison  of 
theory  to  the  experimental  data. 

To  sort  out  the  problem,  we  chose  two  cases,  thick  40nm  well 
(bulk)  and  thin  7.5  nm  QW.  In  thick  well  all  the  complicated 
phenomena  related  to  the  QW  vanish  and  only  CS  is  expected  to  be 
fundamental.  We  investigated  the  CS  effect  on  Jth  using  standard 
bulk  theory  [Ster73].  In  thin  7.5  nm  QW  where  quantum  size  and 
nonradiative  effects  are  important,  we  studied  the  effect  of  CS  on  Jth 
using  valence  band  mixing  (VBM)  QW  theory.  In  the  new  model,  the 
leakage  current  was  included  into  the  total  current  in  terms  of 
spontaneous  emission  efficiencies  Cnsp),  which  can  be  expressed  as 
the  product  of  two  efficiencies,  r\in  = injection  efficiency  and  rir  = 
active  layer  radiative  efficiency  [Zory91,  Lim91].  Injection  efficiency 
was  estimated  from  the  differential  quantum  efficiency  versus  laser 
cavity  length  measurements.  Radiative  recombination  rate  was 
calculated  theoretically.  To  calculate  the  band  structure  and  electron 
(or  hole)  wavefunctions  in  a quantum  well,  we  used  the  multiband 
effective-mass  theory  (k.p  method)  [Lutti55] . For  simplicity,  we 
assumed  that  the  conduction  and  valence  bands  are  decoupled,  which 
is  in  fact  a good  approximation  for  wide-band-gap  semiconductors 
such  as  GaAs  [Sand87].  For  the  conduction  band,  a simple  isotropic 
parabolic  band  was  used.  For  the  valence  band,  we  used  a 4x4 
Luttinger-Kohn  Hamiltonian,  the  elements  of  which  are  derived 
based  on  the  k.p  method.  The  multiband  effective  mass  equations 
were  solved  numerically  using  the  finite-difference  method. 
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From  the  above  work  we  found  that  a reasonable  fit  to  the 
thick  QW(bulk)  experimental  data  at  both  25°C  and  125°C  could  be 
obtained  using  strict  k-selection  bulk  theory  and  a knowledge  of  the 
carrier  injection  efficiency(rnn).  In  making  the  fit,  carrier  scattering 
is  not  necessary,  quite  contrary  to  Yamada's  results  [Yama81].  This 
finding  suggested  the  possibility  that  the  effect  of  CS  on  Jth  might  not 
be  important  even  in  thin  QW  lasers.  This  was  verified  with  VBM  QW 
theory  combined  with  a knowledge  of  the  carrier  injection 
efficiency(rtjn). 

The  thesis  is  organized  as  follows.  We  start  with  a comparison 
of  predicted  Jth  and  T0  with  experimental  data  in  chapter  2.  In 
Chapter  3 we  present  the  results  of  fitting  thick  well  data.  In  Chapter 
4 we  summarize  the  results  of  fitting  thin  QW  data.  Finally,  we  end 
with  summary  and  conclusions  in  chapter  5. 


CHAPTER  2 

COMPARISON  OF  EXPERIMENTAL  J*,  T0  WITH  CALCULATED  Jp,,  T0 

2.1.  Laser  Material 

The  lasers  material  used  in  these  experiments  was  separate 
confinement  heterostructure(SCH)  material  with  a 7.5,  20,  40  nm 
thick,  single-quantum- well(SQW)  GaAs  active  layer  sandwiched 
between  two  AlGaAs  layers  with  A1  composition  graded  between  30% 
and  70%.  The  GRIN-SCH-SQW  material  [Tsan81]  shown  in  Fig.  2-1  was 
grown  in  conventional  metal  organic  chemical  vapour  deposition 
(MOCVD)  vertical  reactor  using  procedures  similar  to  those  described 
by  Dapkus  [Dapk82].  The  specifics  of  the  seven  layer  wafer,  grown  at 
a susceptor  temperature  of  750°C  on  an  n-type  GaAs  substrate(Si 
= 1018/cc),  are  as  follows:  a 1.5-pm  Alo.7Gao.3As  n-type  confining 
layer(Si  = 1018/cc),  a 0.2-pm  AlxGai-xAs  n-type  graded  (x=0.7-0.3) 
barrier  layer(Si  =1017/cc),  an  undoped  GaAs  active  layer,  a 0.2-pm 
AlxGai_xAs  p-type  graded  (x=0. 3-0.7)  barrier  layer(Mg  = 1017/cc),  a 
1.5-pm  Alo.7Gao.3As  p-type  confining  layer  (Mg  =1018/cc),  and  a 0.2- 
pm  GaAs  Mg-doped  GaAs  p+ cap  layer  (Mg  =1019/cc). 

After  growth,  150  pm  conducting  stripes  were  defined  on  the 
material  using  0.15  pm-thick  SiC>2.  For  ohmic  contact,  metal 
deposition  was  made.  The  thickness  of  gold  layer  is  2000  A for  p- 
side  and  1500  A for  n-side.  But  more  gold  was  needed  for  good  wire 
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Fig.  2-1.  GRIN-SCH-SQW  laser  structure. 
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bonding  on  p-side,  and  for  heat  spreader  on  n-side.  So  additional 


deposited.  To  make  sure  of  good  adhesion  of  the  additional  gold  layer 
on  the  existing  metal,  500  A of  Ti  was  deposited  first.  The  deposition 
was  done  with  Electron-Beam  Evaporator  and  the  pressure  during 
the  deposition  was  maintained  at  5xl0‘6  Torr.  Bars  of  various  length 
were  then  cleaved,  separated  into  250,  500,  750,  1000  pm-wide 
chips  and  soldered  p-side  up  on  the  copper  blocks. 

2.2.  Modal  property 

Since  the  stimulated  emission  and  gain  coefficient  are  related 
to  the  photon  flux,  it  is  important  to  consider  the  field  distribution  in 
order  to  determine  the  fraction  of  the  optical  mode  within  the  laser 
active  layer.  The  extent  of  light  coefficient  within  the  active  layer 
may  be  represented  by  the  confinement  factor  T.  Therefore  T is 
defined  as  the  ratio  of  the  light  intensity  within  the  active  layer  to 
the  sum  of  light  intensity  both  within  and  outside  the  active  layer, 
which  can  be  given  by  [Bote84] 


gold  layers  of  2000  A on  p-side  and  4000  A on  n-side  were 


r = 


d/2 

J I Ey | 2 dZ 
-d/2 


(2.1) 


where  Ey  the  electric  field  in  the  direction  perpendicular  to  the 
junction  plane. 


1 0 


GRIN-SCH  waveguide  consists  of  a high-index  low-energy-gap 
active  layer  sandwiched  between  optical  confining  layers,  which  are 
in  turn  surrounded  by  waveguide  cladding  layers.  The  index  profile 
of  the  GRIN-SCH  waveguide  is  shown  in  Fig.  2-2.  The  confining  layers 
generally  have  a nonuniform  index  of  refraction  which  might  vary  in 
a parabolic  or  linear  manner.  The  square  of  the  index  of  refraction  is 
assumed  to  be  linear  in  the  confining  layer.  This  form  is  used  so  that 
Airy  function  solutions  may  be  used  in  the  graded  region.  The 
differences  between  a square  of  the  index  and  a linearly  varying 
index  should  be  negligible  for  most  practical  situations,  since  the 
index  differences  between  layers  are  small. 

Since  most  multiheterostructure  lasers  operate  with  electric 
fields  in  the  plane  of  the  layers,  we  will  calculate  only  such  TE 
modes.  The  general  field  equation  is 

dF,2 

dz2  + (n2k2  - (32)E  = 0 (2.2) 

where  a common  time-  and  space-dependent  factor  of  expfj(cot-px)] 
has  been  removed.  The  x-directional  propagation  coefficient  p is  the 
eigenvalue  to  be  found,  and  the  frequency  co  = 2nc/X.  The  index  in 
the  graded  regions  indicated  in  Fig.  2-2  are 

nz  = (n32  - n22)(z-d/2)/w  + n22  d/2  < z < d/2+w, 

(n22  ' n32)(z+w+d/2)/w  + n32  -w-d/2  < z < -d/2  (2.3) 
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Fig.  2-2.  Spatial  profile  of  the  square  of  the  index  of  refraction  for 
a GRIN-SCH  waveguide.  nl=ni2,  n2=n22,  and  n3=n32 


The  non-normalized  fields  in  the  three  regions  are 


Ey  = Fa  exp(iXaz)  +Ga  exp(-ilaz) 
Ey  = Fb  Ai(t)  + Gb  Ai(t) 

Ey  = Gc  exp[-xc(z-d/2-w)] 


-d/2<  z < d/2, 
d/2<  z < d/2+w 
d/2+w  < z (2.4) 


The  factors  in  the  arguments  of  Eq.(2.4)  are 

Xa  = (nj2k2  - p2),  Xc  = (p2  . ni2k2), 

t=(j)i/3(z-C),  (J)=(n22  - n32)k2/w,  £=  (<})d/2  + n22k2  - p2)/(j)j 


and  the  functions  Ai  and  Bi  are  the  standard  Airy  functions  [Abra72]. 
The  boundary  conditions  are  the  TE  fields  and  their  derivative 
should  be  continuous  at  z = d/2  and  z = w + d/2.  Derivative  of  the 
field  at  z = 0 are  zero.  When  these  boundary  conditions  are  applied  to 
Eq.(2.4),  this  produces  five  equations  linear  in  the  unknown 
coefficients  Fa,  Fb,  Ga,  Gb,  and  Gc,  with  functional  arguments 
dependent  on  the  coefficient  (3.  After  eliminating  Fa,  Fb,  Ga,  Gb,  and 
Gc,  the  resulting  eigenvalue  equation  is 


where  FI  = -Ai’(ac)  Bi(cc)  + Ai(cc)  Bi'(ac), 

F2  = ((j)l/3/Xc)[Ai,(cc)Bi,(ac)  - Ai’(ac)Bi'(cc)], 
F3  = Ai(ac)Bi'(cc)  - Ai'(cc)Bi(ac)], 

F4  = (xc/4>  1/3)[Ai(ac)Bi(cc)  - Ai(cc)Bi(ac)], 
ac=  (})1/3(d/2-  Q,  and  cc  = (|)  1/3(d/2+w-0 


k = 2t i/X, 


(2.5) 


(2.6) 


(2.7) 


The  primed  functions  are  spatial  derivatives.  In  numerical 
calculation,  we  used  IMSL  libraries  AI,  BI,  AID,  and  BID  for  the  Airy 
functions  and  its  derivatives.  The  parameters  used  for  this  structures 
are  ni=3.62)  n2=3.4  , n3=3.1  [Case78],  w=0.2  pm,  and  X=0.88  pm.  In  Fig. 
2-3  is  shown  plot  of  confinement  factor  T as  a function  of  active 
layer  thicknesses  from  5 nm  to  40  nm. 

2.3.  Measurement 

Three  different  wafers  with  the  structure  shown  in  Fig.  2-1 
were  used.  The  only  variable  between  the  three  wafers  was  the 
active  layer  thickness  d which  was  7.5  nm  (SQW  lasers)  , 20  nm,  and 
40  nm  (bulk  lasers).  Oxide  stripe  lasers  were  fabricated  from  each 
wafer  with  a stripe  width  W of  150  pm  and  cavity  length  L values 
ranging  from  250  pm  to  1000  pm.  The  facets  were  uncoated  and  only 
lasers  with  uniform  near  fields  were  used.  The  measurement  set-up 
is  shown  schematically  in  Fig.  2-4.  The  samples  were  measured  with 
pulsed  power  supply.  With  the  copper  block  held  at  a certain 
temperature,  output  power/current  profiles  were  generated  using 
200  nm-wide  drive  pulses.  Temperature  was  controlled  by  changing 
the  voltage  applied  to  the  heating  element  and  monitored  by 
temperature  indicator. 

For  a precise  determination  of  threshold  current  (Ith)>  the 
spectrum  of  each  laser  was  observed  as  a function  of  current  (low 
duty  cycle  pulsed  operation).  The  differential  quantum  efficiency  (qd) 
of  each  laser  was  measured  from  the  slope  of  the  output  power  (P0) 
vs.  current  (I)  characteristic.  Internal  loss  coefficient  aj  and  carrier 

i 

injection  efficiency  r|jn  shown  in  Table  I were  determined  from  rid  vs. 


Confinement  factor 
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QW  thickness  (nm) 


Fig.  2-3.  Optical  confinement  factor  r as  a function  of  active  layer 
thickness  d for  the  GRIN-SCH  trangular  well  waveguide. 
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Fig.  2-4.  The  schematic  of  the  threshold  current 


measurement  set-up 
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L measurements  [Zmud91].  The  measured  values  of  Jth  are  given  in 
Table  II  [Zmud91].  Figs.  2-5  and  2-6  show  the  rj 1 vs  L measurement 
data  for  the  40  nm  and  7.5  nm  QW  respectively.  The  peak  of  the 
mode  gain  at  threshold  is  given  by  the  total  cavity  loss  Gth  = rgth  = cq 
+ j~ln~  , where  R is  the  mirror  reflectivity  (0.32)  at  facet,  and  L is  the 

device  length.  The  material  peak  gain  at  threshold,  gth,  was 
calculated  using  the  relation  between  mode  gain  and  active  layer 
material  gain,  Gth  = Ggth- 


2.4.  Prediction  of  J^and  T n 

For  the  prediction  of  the  threshold  current  density  and 
characteristic  temperature  at  25°C,  a detailed,  published  CZR  model 
[Chin88,  Chin89]  was  used.  Important  features  included  in  the  Chinn 
model  are  as  follows:  carrier  scattering,  gain  anisotropy,  V band  non- 
parabolicity,  energy  levels  of  GRIN  layer,  L and  X type  carriers,  band 
gap  renormalization,  carrier  leakage  from  GRIN  layer,  light  hole 
subband  nonparabolicity,  nonradiative  processes  by  carrier  leakage 
over  the  confining  potential  barriers  and  from  the  Auger 
recombination.  For  the  prediction  of  threshold  current  density,  we 
have  used  the  Gth  estimated  in  Table  II  for  a given  d and  L. 

The  calculation  of  the  peak  gain  requires  a calculation  of  the 
spectral  dependence  of  the  mode  gain  from  the  gain  equation, 
followed  by  spectral  convolution  to  include  the  intraband  carrier 
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Table.  I 

Measured  values  of  r|  in,  aj,  and  calculated  values  of 
T for  lasers  from  each  wafer 


d (nm)  T|in 

cci(cnr 

!)  r 

7.5 

0.82 

7.0 

0.030 

20 

0.90 

8.5 

0.081 

40 

0.98 

13 

0.180 

Table.  II 

Measured  values  of  Jth  and 

calculated  values  of  gth  and  Jth/d 

(nm) 

L(fim) 

Jth  (A/cm2) 

gth(cm-l)  Jth/d  (kA/cm2-|im)  Gth  (cm-1) 

7.5 

250 

448 

1752 

59.7 

52.6 

if 

500 

239 

992 

31.9 

29.8 

it 

750 

178 

738 

23.7 

22.2 

it 

1000 

161 

612 

21.5 

18.4 

20 

250 

452 

654 

22.6 

54.1 

II 

500 

323 

378 

16.2 

31.3 

it 

750 

245 

286 

12.3 

23.7 

ii 

1000 

212 

240 

10.6 

19.9 

40 

250 

467 

335 

11.7 

60.6 

it 

500 

371 

205 

9.28 

37.1 

it 

750 

319 

161 

7.98 

29.1 

ii 

1000 

277 

139 

6.93 

25.2 

L ((jm  ) 


Fig.  2-5.  Reciprocal  of  the  external  differential  quantum  efficiency 
(1/rid)  vs  cavity  length  for  40  nm  GRIN-SCH-SQW  at  25°C. 


L ((jm) 


Fig.  2-6.  Reciprocal  of  the  external  differential  quantum  efficiency 
(1/rid)  vs  cavity  length  for  7.5  nm  GRIN-SCH-SQW  at  25°C. 


20 


scattering.  The  sharp-featured  curve  represents  the  unconvolved 
spectra,  and  the  corresponding  curve  of  the  dotted  line  type 
represents  the  convolved  spectra  in  Fig  2-7.  The  model  gain  has  been 
calculated  using  the  product  of  the  optical  confinement  factor  times 
the  material  gain  coefficient  evaluated  at  the  curve's  spectral  peak. 
Each  curve  is  associated  with  a well-defined  carrier  density. 

In  order  to  make  a correspondence  between  carrier  density 
and  measurable  current  density  (proportional  to  carrier  generation 
rate),  we  must  include  all  of  the  subthreshold  carrier  recombination 
processes,  sum  their  rates,  and  equate  them  to  the  carrier  injection 
rate.  The  radiative  component  of  the  excitation  current  density  at  the 
various  drive  levels  is  found  from  the  integral  over  energy  of  the 
total  spontaneous  emission  spectrum.  An  example  of  such  a 
calculated  spectrum  is  shown  in  Fig. 2-8,  broadened  by  the  addition 
of  intraband  scattering.  Nonradiative  components  of  the  current 
density  include  contributions  from  indirect-gap  nonradiative 
recombination,  Auger  recombination,  and  leakage  of  carriers  over  the 
confining  energy  barriers.  The  radiative  and  nonradiative 
recombination  rates  can  all  be  calculated  for  each  value  of  carrier 
density.  Since  the  peak  gain  is  a function  of  carrier  density,  each 
value  of  peak  gain  can  be  associated  with  a current  density  whose 
components  can  be  added  sequentially  for  modeling  purposes,  to 
examine  their  relative  sizes.  The  carrier  density,  while  a necessary 
variable  for  the  model,  does  not  appear  in  the  final  result  for  the 
peak  mode  gain  G as  a function  of  current  density  J (Fig.  2-9).  T0  can 
be  predicted  simply  by  calculating  G(J)  at  several  temperatures,  and 


Energy  (eV) 


Fig.  2-7.  TE  mode  gain  as  a function  of  energy  for  a GRIN-SCH-SQW 
laser.  Unconvolved:  dashed  line  and  convolved:  dotted  line. 
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Fig.  2-8.  Calculated  total  spontaneous  emission  spectral  density 
(1015  cm'2  sec_1  nm*1)  as  a function  of  wavelength.  The 
spectrum  has  been  convolved. 
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Fig.  2-9.  Peak  TE  mode  gain  as  a function  of  current  density.  The 
dashed  curves,  from  unconvolved  spectral  gains,  are 
calculated  for  radiative  recombination(R),  radiative  plus 
Auger  recombination(R-t-A),  and  radiative  and  Auger 
recombination  plus  carrier  leakage(R+A+L)  over  the 
potential  barriers.  The  corresponding  solid  curves  come 
from  spectrally  convolved  gains.  CS  is  carrier  scattering. 


finding  the  required  threshold  current  density  at  those  temperatures 
for  a specified  threshold  gain.  The  T0  was  obtained  from  the 
exponential  variation  of  Jth  with  temperature,  J(T)=J(Tr)exp[(T- 
Tr)/T0L  where  Tr  is  the  arbitrary  reference  temperature,  in  the 
range  25°C  - 70°C.  The  predicted  peak  mode  gain  vs  current  density 
and  T0  vs  QW  thickness  are  shown  in  Figs.  2-10  and  2-11. 

2.5.  Results  and  discussion 

We  compared  the  measured  Gth  vs  Jth  with  predicted  one  for  the 
d=7.5  nm,  20  nm  and  40  nm  QW  in  Fig. 2-12.  In  20  nm  and  40  nm 
QW,  theoretical  Jth  is  uniformally  overestimated,  while  in  7.5  nm, 
experimental  values  are  close  to  the  prediction  for  long  lasers  (750 
and  1000  pm).  However,  theory  starts  to  deviate  from  the 
experiment  as  L shortens  (250  and  500  pm).  Threshold  current 
density  as  a function  d for  L=500  pm  are  shown  in  Fig. 2-13.  For  the 
theoretical  mode  gain  estimation,  we  used  aj  values  in  the 
intermediate  well  size  interpolated  linearly  from  the  measured  aj 
values.  In  Fig. 2-14  we  plotted  characteristic  temperature  (T0)  as  a 
function  of  quantum  well  thickness  d for  these  structures.  Predicted 
T0  values  agree  pretty  well  with  experimental  values  but  are 
overestimated  because  Jth  is  not  predicted  properly.  There  is  a 
general  tendency  for  T0  to  increase,  as  shown  in  Fig.  2-14,  as 
quantum  well  thickness  increases  [Hers84].  Thick  QW  diode  lasers 
are  favored  for  the  high  temperature  operation  from  the  T0  point  of 
view.  T0  is  a very  complicated  parameter  as  found  by  several 
authors.  Experimental  work  indicates  that  T0  is  actually  decreasing 
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Fig.  2-10.  Peak  TE  mode  gain  as  a function  of  current  density  for  a 
GRIN-SCH-SQW  laser  as  a function  of  quantum  well 
thickness. 
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Fig.  2-11.  Characteristic  temperature  of  GRIN-SCH-SQW  laser  as 
function  of  quantum  well  thickness. 
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Fig.  2-12.  Comparison  of  measured  and  theoretical  mode  gain  vs 

current  density  for  a GRIN-SCH-SQW  laser  as  a function  of 
quantum  well  thickness  . 
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Fig.  2-13.  Comparison  of  measured  and  theoretical  mode  gain  vs 
current  density  as  a function  of  QW  thickness  for  GRIN- 
SCH-SQW  laser  with  cavity  length  L=500  pm  and  stripe 
width  W=150  pm  diode. 
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Fig.  2-14.  Comparison  of  measured  and  theoretical  characteristic 

temperature  (T0)  as  a function  of  QW  thickness  for  cavity 
length  L=500  pm  and  stripe  width  W=150  pm  diode  laser. 
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function  of  temperature  [Zory86a,  Shea87].  Also  it  has  been  found 
that  T0  for  a GRIN-SCH-SQW  increases  with  cavity  length  [Zory86b, 
Leop87],  T0  is  very  bad  parameter.  Therefore  we  do  not  study  T0 
anymore  in  this  thesis. 

The  possible  reasons  why  the  CZR  model  do  not  agree  with 
experiment  are  as  follows:  (i)  CZR  model  breaks  down  for  the  thick 
well  (bulk),  (ii)  Lineshape  function,  which  was  assumed  to  be 
independent  of  well  thickness,  temperature,  carrier  energy  and 
carrier  density,  is  not  appropriate,  (iii)  Nonradiative  current  might 
not  be  estimated  correctly,  (iv)  CZR  model  does  not  treat  the  valence 
band  correctly.  Fundamental  studies  [Cola87,  Ahn88]  have  shown 
that  the  valence  band  structure  of  quantum  wells  deviate 
significantly  from  that  of  the  parabolic  band  due  to  VBM  as  the  QW 
width  decreases.  VBM  gives  nonparabolic  bands  which  affect  both 
the  minimun  confinement  energies  and  the  density  of  states 
functions.  The  density  of  states  functions  are  found  to  differ 
considerably  from  the  the  simple  step-like  shape  computed  from  the 
band  edge  effective  masses.  The  band  mixing  also  results  in  large 
differences  in  the  momentum  vector  dependence  of  the  matrix 
elements,  and  moreover,  some  of  the  transitions  which  were 
previously  assumed  to  be  forbidden  become  partially  allowed.  The 
quantum  well  gain  calculated  by  the  k.p  method  and  k-selection 
rules  clearly  show  the  effects  of  band  mixing  both  in  shape  and  in 
peak  magnitude.  These  results  are  considerably  different  than  those 
obtained  from  the  conventional  methods  with  parabolic  bands  and 
with  matrix  elements  which  have  constant  total  magnitude.  The 
threshold  current  in  a QW  laser  is  found  to  be  strongly 


underestimated  if  band  mixing  effects  are  ignored.  Although  the 
effects  of  the  band-mixing  on  the  band  structure  and  gain  of  a QW 
have  been  studied  [Cola87,  Ahn88,  Zory92],  very  little  work  has  been 
done  on  the  comparison  of  theory  with  the  experimental  data. 

2.6.  Suggestion 

To  sort  out  the  problem,  we  chose  two  cases,  thick  40nm  well 
(bulk)  and  thin  7.5  nm  QW.  In  thick  well  all  the  complicated 
phenomena  related  to  the  QW  vanish  and  only  CS  is  expected  to  be 
fundamental.  We  investigated  the  CS  effect  on  Jth  using  standard 
bulk  theory  [Ster73].  In  thin  7.5  nm  QW  where  quantum  size  and 
nonradiative  effects  are  important,  we  studied  the  CS  effects  on  Jth 
using  VBM  QW  theory.  In  the  new  model,  the  leakage  current  was 
included  into  the  total  current  in  terms  of  spontaneous  emission 
efficiencies  (r|Sp),  which  can  be  expressed  as  the  product  of  two 
efficiencies,  rjin  = injection  efficiency  and  r|r  = active  layer  radiative 
efficiency  [Zory91,  Lim91]. 


CHAPTER  3 

FITTING  THICK  QW  DATA  WITH  STANDARD  BULK  THEORY 

3.1.  Introduction 

Threshold  current  density  Jth  predictions  for  diode  lasers  are 
based  on  a knowledge  of  the  relationship  between  mode  gain  G and 
current  density  J for  the  laser  of  interest.  In  order  to  demonstrate 
the  validity  of  a given  G-J  relation,  one  usually  performs  experiments 
on  simple,  wide  stripe  lasers  where  complications  due  to  current 
spreading,  lateral  mode  overlap  and  antiguiding  can  be  minimized.  In 
this  case,  G can  be  expressed  as  the  product  of  a transverse  mode 
confinement  factor  T and  peak  material  gain  g.  Mode  gain  at 
threshold,  Gth,  can  then  be  expressed  in  the  form  [Chin88], 

Gth  = r gth  = ai  + LlnIC  (3.1) 

where  aj  is  the  mode  loss  coefficient,  L is  the  cavity  length  and  R is 
mode  facet  reflectivity.  Since  one  can  estimate  T,  a*  and  R with  fair 
accuracy  for  a given  resonator,  one  can  estimate  gth. 

To  find  the  Jth  value  associated  with  the  estimated  gth  value, 
one  usually  uses  a strict  k selection  theory  and  a carrier  scattering 
(CS)  convolution  procedure  [Chin88].  If  one  knows  the  CS  function,  one 
can  calculate  the  value  of  carrier  density  N which  yields  gth.  This 
value  of  N = Nth  is  then  used  to  compute  the  spontaneous 
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recombination  rate  density  at  threshold  (Rsp)th  which  is  related  to  Jth 
by 


Jth 


ed(Rsp)th 

Tlsp 


(3.2) 


In  this  equation,  e is  the  electron  charge,  d is  the  active  layer 
thickness  and  qsp  is  the  efficiency  for  converting  J into  spontaneous 
emission  in  the  active  layer. 

In  experimental  situations  where  threshold  gain  requirements 
are  high  (for  example  when  the  active  layer  temperature  T is  high,  L 
is  short  or  d is  very  thin),  it  is  difficult  to  estimate  r|Sp  [Chin88]. 
However,  if  we  choose  T low,  L long  and  d thick,  then  rj  sp  should  be  = 
1,  provided  the  laser  material  being  used  is  good.  The  MOCVD  laser 
material  [Dapk82]  used  in  this  work  is  shown  in  Fig.  2-1.  It  is  a linearly 
graded  refractive  index  (GRIN),  separate  confinement  heterojunction 
(SCH),  single  potential  well  (SW)  structure  with  an  active  layer 
thickness  d = 40  nm.  The  GRIN  layers  on  either  side  of  d are  0.2  pm 
thick  with  aluminum  grading  from  30%  to  70%.  The  p and  n doped 
cladding  layers  have  70%  aluminum  and  are  1.5  pm  thick.  The  value 
of  G for  this  structure  is  the  product  of  the  geometrical  overlap  factor 
(0.164)  times  a correction  factor  of  1.1. 


3.2.  Predictions 

Fig.  3-1  and  3-2  show  calculated  plots  of  gth  vs  Jth/d  at  T = 25 
°C  and  125  °C,  respectively,  without  CS  convolution  ( = g0)  and  with  CS 
convolution  ( = gc).  The  convolution  function  (L)  used  in  this  case  is 


Threshold  Gain  gth  (cm 
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Jth/d  (kA/cm2  -jim) 


Fig.  3-1.  Threshold  gain  as  a function  of  Jth/d  at  25°C.  Solid  curve  is 
without  CS  convolution  (g0)  and  dashed  curves  is  with  CS 
convolution  (gc).  m*  in  matrix  element  is  0.053mo.  T|sp  is 
assumed  to  be  1. 
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Jth/d  (kA/cm2’  -|im) 


Fig.  3-2.  Threshold  gain  as  a function  of  Jth/d  at  125°C.  Solid  curve 
is  without  CS  convolution  (g0)  and  dashed  curves  is  with  CS 
convolution  (gc).  m*  in  matrix  element  is  0.053mo.  rjsp  is 
assumed  to  be  1. 
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Lorentzian  with  a width  Ycv  which  depends  on  photon  energy  hv  and 
active  layer  temperature  T [Yama81,  Asad89]  : 


L = L(hv 


Ycv 

(hv-Ecv)2  + y2cv 


(3.3) 


where  Ycv  = Ycv(hv,T)  = ~ — . Band  edge  scattering  times  xcv(hv  = 

tcv  (hv, 1 ) 

Eg,  T)  of  0.1  ps  at  25  °C  and  0.055  ps  at  125  °C  were  used  to  compute 
the  gc  curves  in  Figs.  3-1  and  3-2.  Ecv  is  the  photon  energy  at  which 
the  convolved  gain  is^  calculated  when  integrating  the  convolution 
integral  over  values  of  hv  from  Eg  to  °°. 

The  expression  for  spectral  gain  g(hv,T)  that  is  convolved  with 
Eq.  (3.3)  to  determine  the  gc  curve  (no  convolution  for  the  g0  curve) 
is  [Yan90]  : 

1 e2h 

*(hv’T)  = IM°|2pri  ^ (3-4) 


where  n is  a refractive  index  parameter,  IM0I2  is  the  transition 
matrix  element  given  by 

IMqI2  (3.5) 

6(Eg+y) 

where  Eg  = energy  band  gap,  A = split-off  valence  band  separation, 
and  prd  is  the  reduced  density  of  states  with  spin  degeneracy  factor 
of  2 included  given  by 
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Prd  — 


8/171111  r3/2 


(hv-Eg)l/2 


(3.6) 


The  value  for  (Rsp)th  used  to  determine  Jth  in  Eq.(3.2)  is  calculated  by 
using  Nth  in  the  following  [Yan90]  : 


oo  oo 

r\ 

(Rsp)th  = J[rSp(hv,T)]th  d(hv)  = -_4^e  n 1^2  f hv  prd  fc(l-fv)d(hv)  (3.7) 
Eg  e0h2c3m02  Eg 


Effective  masses  used  in  these  calculations  [Corz90]  were  mc  = 
0.067mo,  mhh  = 0.59mo  and  mjh  = 0.083mo.  The  value  of  m*  used  in 
Eq.  (3.5)  was  0.053mo  [Yan90]  and  the  value  of  n and  n'  in  Eqs.  (3.4) 
and  (3.7)  was  3.6.  Split-off  valence  band  separation  was  0.33  eV. 

It  is  interesting  to  note  that  if  we  change  T|sp  in  the  25  °C  g0 
curve  of  Fig.  3-1  from  1 to  0.61,  it  closely  approximates  the 
convolved  curve  (gc)  where  r|sp  = 1.  For  the  125  °C  curve,  qsp  = 0.54 
achieves  the  same  result.  It  is  clear,  therefore,  that  one  must  know 
lisp  if  one  is  to  utilize  this  method  to  determine  the  correct  CS 
convolution  function.  As  mentioned  earlier,  we  chose  to  use 
relatively  thick  active  layer  material  in  this  work  (d  = 40  nm)  so  that 
we  could  make  the  r|sp  = 1 assumption  with  some  confidence  when 
comparing  measured  Jth  data  of  long  lasers  at  low  temperature  with 
theory.  The  validity  of  this  assumption  is  discussed  below. 
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3.3.  Results  and  discussions 

In  the  previous  section,  we  showed  that  the  G-J  relation 
when  applied  to  the  wide  stripe  lasers  described  earlier,  can  be 
transformed  to  the  gth  - Jth/d  relations  shown  in  Figs.  3-1  and  3-2. 
Fig.  3-3  is  the  same  plot  as  Fig.  3-1  with  the  experimental  data  points 
added.  The  experimental  gth  values  shown  are  determined  using  Eq. 
(3.1)  and  measured  aj  values  deduced  from  differential  quantum 
efficiency  measurements  vs  L.  The  error  bars  on  the  gth  points 
represent  the  possible  variation  in  gth  due  to  uncertainties  associated 
with  G,  aj  and  R.  The  error  bars  on  the  Jth/d  points  represent  the 
possible  variation  in  Ith  from  chip  to  chip  and  uncertainty  in  amount 
of  current  spreading.  The  lowest  point  represents  the  L = 1000  pm 
lasers  and  the  highest  represents  the  L = 250  pm  lasers.  As  shown, 
the  calculated  g0  curve  overestimates  the  value  of  Jth  by  about  8% 
whereas  the  gc  curve  overestimates  Jth  by  about  78%.  In  this  case,  it 
is  impossible  to  fit  the  g0  curve  to  the  experimental  points  by  making 
T|sp  < 1,  since  this  approach  moves  the  g0  curve  to  the  right.  It  is  also 
impossible  to  fit  the  gc  curve  to  the  data  since  gc  is  always  to  the 
right  of  g0.  If  the  expression  for  IM0I2  in  Eq.(3.5)  is  modified  by  using 
m*  = 0.067mo  rather  than  0.053mo,  the  recalculated  g0  curve  (g0') 
shown  in  Fig.  3-4  is  to  the  left  of  the  experimental  points.  In  this 
case,  if  one  uses  rj sp  = 0.98  in  the  expression  for  Jth  (Eq.(3.2)),  g0'  fits 
well  with  the  data.  The  gc  curve  of  course  has  no  chance  to  fit  the 
data  unless  tisp  > 1,  prohibited  by  definition.  In  order  to  determine  if 
T|sp  = 0.98  is  reasonable,  we  utilized  the  fact  that  risp  can  be 
expressed  as  the  product  of  two  efficiencies,  ri  jn  = injection  efficiency 
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Jth/d  (kA/cmz  -jim) 


Fig.  3-3.  Threshold  gain  as  a function  of  Jth/d  at  25°C  with 

m*=0.053mo  in  matrix  element.  Solid  curve  is  without  CS 
convolution  (g0)  and  dashed  curves  is  with  CS  convolution 
(gc).  "Hsp  is  assumed  to  be  1. 


Threshold  Gain  gth  (cm"1) 
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Jth/d  (kA/cmz  -jim) 


Fig.  3-4.  Threshold  gain  as  a function  of  Jth/d  at  25°C  with 

m*=0.067mo  in  matrix  element.  Solid  curve  is  without  CS 
convolution  (g0)  and  dashed  curves  is  with  CS  convolution 
(gc).  Tlsp  is  assumed  to  be  1. 
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and  r|r=  active  layer  radiative  efficiency  [Zory91].  From  differential 
quantum  efficiency  (rjd)  vs  L measurements,  we  determined  that  Hin 

riSpis  assumed  to  be  1=  0.98,  implying  that  tl  r = 1.  Since  r\  r = 

Rsp/(Rsp+Rnr)  where  Rnr  is  the  non-radiative  recombination  rate 
density  in  the  active  layer  at  threshold,  T|r=  1 implies  Rnr  = 0,  not 

unreasonable  for  high  quality  diode  laser  material  at  low  carrier 
densities. 

The  justification  for  making  the  particular  modification  to  IM0I2 
described  above  is  essentially  historical  in  that  the  resulting  g0' 
curve  is  the  same  as  would  be  predicted  using  Stern's  k-selection 

formalism  [Ster73].  In  reviewing  how  well  Stern's  theory  compared 

to  experiment,  Thompson  [Thom80]  pointed  out  that  it  predicted  Jth 
values  somewhat  lower  than  that  measured  for  their  best  lasers. 
About  the  same  time  (early  ’80s),  Yamada  pointed  out  that  Stern's 
theory  did  not  take  CS  into  account  implying  that  Stern's  g(hv) 
function  was  too  optimistic.  By  convolving  the  g(hv)  function  with  a 
Lorentzian  broadening  function  (CS  time  of  about  10' 13  sec),  Yamada 
and  Ishiguro  were  able  to  fit  their  data  [Yama81].  On  the  other  hand, 
Garbuzov  et  al  [Garb87]  showed  good  agreement  with  Jth  values 
measured  by  Hakki  and  Paoli  with  a carrier  scattering  time  of  about 
10’ 12  sec  (CS  contribution  small).  The  essence  of  the  above  history 
appears  to  be  that  if  measured  Jth  values  are  low  enough,  Stern's 
theory  without  CS  is  sufficient  to  explain  results. 

Fig.  3-5  shows  the  calculated  g0'  curve  at  T = 125  °C  for  r|sp  = 1 
along  with  the  measured  125°C  data.  Since  the  data  points  are  to  the 
right  of  the  predicted  curve,  one  can  fit  the  data  by  using  CS,  tisp  < 1 
or  a combination  of  both.  Fig.  3-6  shows  that  if  we  choose  Tisp  = 0.74, 
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Jth/d  (kA/cm^-jim) 


Fig.  3-5.  Threshold  gain  as  a function  of  Jth/d  without  CS 

convolution  at  125°C  with  risp  =1.  m*  in  matrix  element  is 
0.067mo. 
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Jth/d  (kA/cmz  -|im) 


Fig.  3-6.  Threshold  gain  as  a function  of  Jth/d  without  CS 

convolution  at  125°C  with  r|sp  =0.74.  m*  in  matrix  element 
is  0.067mo. 
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the  L = 1000  pm  point  falls  on  the  g0'  curve.  From  the  slope 
efficiency  data  for  lasers  with  L = 500,  750  and  1000  pm  at  T = 25  °C, 
we  determined  that  T|in  = 0.78  at  125°C  (see  Appendix  A),  very  close 
to  the  value  for  qsp  (=0.74)  required  to  fit  the  L = 1000  pm  data 
point.  In  this  case,  g0'  "predicts"  the  measured  Jth  values  using  rj in  = 
0.78,  rjr  = 0.95  for  the  1000  pm  laser  and  decreasing  T|r  for  the 
shorter  lasers.  This  dependence  of  T|r  on  gth  is  very  reasonable  since 
nonradiative  recombination  (e.g.  Auger  recombination)  can  be 
expected  to  lower  rir  as  gth  and  Nth  increase.  Fig.  3-7  shows  the  g vs 
Jth/d  with  fixed  r|jn  = 0.78  and  variable  rir.  Radiative  efficiency  was 
calculated  from  the  ratio  of  the  radiative  recombination  rate  and 
Auger  recombination  rate.  A better  fit  to  the  data  may  be  obtained  if 
exact  riin  for  each  cavity  length  lasers  are  used.  Based  on  these 
results,  it  appears  that  if  facet  coatings  had  been  used  to  lower  the 
required  Nth  for  these  lasers,  then  it  would  have  been  a good 
approximation  to  use  r|r  = 1 for  the  500  and  750  pm  long  lasers,  even 

at  125°C.  In  that  case,  we  would  then  only  have  to  know  how  to 
predict  the  value  of  rj jn  at  125°C  in  order  to  predict  Jth  at  125  °C. 
Since  rj in  is  related  in  a simple  way  to  rid  for  lasers  where  output 
power  vs  drive  current  is  a straight  line  [Zory91],  the  ability  to 
predict  T^n  also  allows  a prediction  of  rid-  Consequently,  if  CS  can  be 
ignored,  then  the  problem  of  predicting  Jth  and  rid  at  any 
temperature  reduces  to  the  problem  of  predicting  the  temperature 
dependence  of  ri in- 

As  mentioned  above,  since  the  125  °C  data  points  of  Fig.  3-5 
are  to  the  right  of  the  g0'  curve  with  iq  sp  = 1,  it  may  be  possible  to 
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J /d  (kA/cm2-um) 


Fig.  3-7.  Threshold  gain  as  a function  of  Jth/d  without  CS 

convolution  at  125°C  with  r\in  =0.78  and  r|r  variable,  m*  in 
matrix  element  is  0.067mo. 
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include  CS  and  still  fit  the  data  since  CS  will  shift  the  curve  to  the 
right.  However,  since  the  measured  value  of  r\ jn  is  0.78,  one  must  set 
T|sp  < 0.78,  which  means  the  difference  between  the  data  points  and 
the  curve  of  Fig.  3-6  actually  indicates  whether  or  not  CS  can  be  used 
to  fit  the  data,  not  the  difference  between  the  data  points  and  the 
curve  in  Fig.  3-5.  Since  including  CS  will  shift  the  curve  in  Fig.  3-6 
substantially  to  the  right,  no  fit  can  be  obtained  at  125  °C  if  CS  is 
included.  Therefore,  it  is  not  possible  to  fit  either  the  25  °C  data  or 
the  125  °C  data  if  CS  is  included.  The  only  way  to  include  CS  and  fit 
the  data  with  allowed  values  of  r)sp  is  to  modify  the  expressions  for 
gain  and  spontaneous  emission  in  Equations  (3.4)  and  (3.7).  For 
example,  if  a factor  of  1/2  is  introduced  on  the  right-side  of  both 
Equations  (3.4)  and  (3.7),  a good  fit  using  the  "accepted"  bandedge  CS 
time  of  0.1  ps  can  be  obtained  at  T = 25  °C  with  qsp  = 1. 

3.5.  Summary 

In  summary,  threshold  current  density  predictions  were  made  for 
GaAs/AlGaAs  GRIN-SCH-SW  diode  lasers  with  40  nm  thick  active 
layers  at  25  °C  and  125  °C.  A reasonable  fit  to  the  experimental  data 
at  both  temperatures  was  obtained  without  including  CS.  Standard 
bulk  theory  combined  with  measured  carrier  injection  efficiency 
(t) in)  was  used. 


APPENDIX  A 

CARRIER  INJECTION  EFFICIENCY  AT  125°C 


Carrier  injection  efficiency  (rjin)  is  determined  from  the  differential 
quantum  efficiency  (rj  d)  vs  L measurements.  However,  it  is  difficult 
to  estimate  hin  in  experimental  situations  where  threshold  gain 
requirements  are  high  (for  example  when  the  active  layer 
temperature  T is  high).  Therefore  we  deduce  hin  from  the  following 
simple  relationship,  r|d  = Tiinr|o(ai),  where  rj 0 is  the  external  quantum 
efficiency  defined  by 


r|0  = 


2 


ajL  + In  — 


(A.l) 


The  differential  quantum  efficiency  ratio  between  125°C  and  25°C 
can  be  rewritten  as  follow  for  facet  mirror  reflectivity  R=0.32 

r|d(1250C)  _ t| in(  1 25°C)  cti(250C)L+1.2 
Tid(25°C)  " r|in(250C)  ai(125°C)L+1.2  (A'2) 


We  assume  that  aj  is  same  for  both  at  25°C  and  125°C.  Then  r|jn  at 
125°C  can  be  represented  as  the  ratio  of  the  diffenertial  quantum 
efficiency  for  both  temperature. 
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Tlin(125°C) 


Tld(125°C) 

rid(250C) 


(A.3) 


Fig.A-1  shows  the  rj  in  at  125°C  as  function  of  cavity  length  L. 


Injection  efficiency 
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Cavity  length  L (pm) 


Fig.  A-l.  Injection  efficiency  as  a function  of  cavity  length  L at  125°C. 


CHAPTER  4 

FITTING  THIN  QW  EXPERIMENTAL  DATA 
WITH  NEW  VBM  QW  THEORY 


4.1.  Introduction 

Fundamental  studies  [Cola87,  Ahn88]  have  shown  that  the 
valence  band  structure  of  quantum  wells  deviates  significantly  from 
that  of  the  parabolic  band  due  to  valence  band  mixing  as  the  QW 
width  decreases.  The  valence  band  mixing  gives  nonparabolic  bands 
which  affect  both  the  subbands  energies  and  the  density  of  states 
functions,  (i)  The  density  of  states  functions  are  found  to  differ 
considerably  from  the  simple  step-like  shape  computed  from  the 
band  edge  effective  masses,  (ii)  The  coupling  of  heavy  and  light  hole 
results  in  large  difference  in  the  matrix  element,  and  some  of  the 
transitions  which  were  previously  assumed  to  be  forbidden  are 
partially  allowed.  This  phenomena  not  predicted  by  the  parabolic 
model  violate  the  An  =0  selection  rule,  where  n is  the  principle 
quantum  number. 

The  quantum  well  gain  calculated  by  the  k.p  method  and  k- 
selection  rules  clearly  show  the  effects  of  band  mixing.  These  results 
are  considerably  different  than  those  obtained  from  the  conventional 
methods  with  parabolic  bands  and  with  matrix  elements  which  have 
constant  magnitude.  The  threshold  current  in  a QW  laser  is  found  to 
be  strongly  underestimated  if  band  mixing  effects  are  ignored. 
Although  the  modification  to  spectral  gain  due  to  band-mixing  has 
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been  discussed  in  the  literature  [Cola87,  Ahn88,  Corz91],  very  little 
work  has  been  done  on  the  comparison  of  theory  with  experiment. 

Using  the  same  procedures  as  used  in  fitting  bulk  data,  we  fitted  the 
thin  7.5  nm  QW  data  with  valence  band-mixing  QW  theory. 

The  multiband  effective-mass  theory  (k.p  method)  [Lutti55] 
was  used  to  calculate  the  band  structure  and  electron  (or  hole)  wave 
functions  of  a quantum  well.  For  simpilicity,  we  assume  that  the 
conduction  and  valence  bands  are  decoupled,  which  is  in  fact  a good 
approximation  for  wide-band-gap  semiconductors  such  as  GaAs 
[Sand87],  For  the  conduction  band,  a simple  isotropic  parabolic  band 
is  used.  For  the  valence  band,  we  used  4x4  Luttinger-Kohn 

Hamiltonian,  the  element  of  which  are  derived  based  on  the  k.p 
method.  We  used  finite  difference  method  to  solve  the  multiband 
effective  mass  equations. 

4.2.  Valence  band  mixing  quantum  well  theory 
4.2.1.  Effective  mass  equations 

The  theory  of  the  k.p  method  in  bulk  semiconductor  was 
introduced  by  Kane  [Kane55]  and  has  been  applied  successfully  to 
bulk  semiconductors  by  number  of  authors  [Lutti55,  Lutt56, 
Lawa71].  The  formalism  was  extended  to  the  study  of  thin  films  by 

S.S.  Nedorezov  [Nedo71],  and  recently  has  been  used  in  the  study  of 

quantum  well  structures  [Faso84,  Sand85,  Chan85,  Cola87,  Ahn88, 
Corz90].  This  approach  has  the  virtue  of  being  relatively  simple  and 
yet  takes  into  account  the  important  aspects  of  the  strong  hole 
interaction  near  the  zone  center.  The  formalism  allows  determination 
of  the  shape  of  the  energy  bands  in  the  neighborhood  of  k=0  through 
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symmetry  considerations  and  an  empirical  fit  of  the  effective  masses. 
In  this  method,  the  wavefunctions  of  the  QW  are  expanded  in  the 
wave  functions  of  both  constituent  materials  around  the  k=0  point. 
The  resulting  Hamiltonian  for  the  slowly  varying  part  of  the  wave 
function,  the  envelope  function,  can  be  described  by  the 
corresponding  k.p  matrix  in  each  constituent  material. 

The  subband  structure  of  a quantum  well  is  obtained  by 
solving  the  Schroedinger  equation. 

p2 

Ho'F  = [ 2^  + V(r)  ] W = EXF’  (4.1) 

where  p is  the  momentum  operator,  H0  is  the  perfect  crystal 
Hamiltonian  and  V(r)  is  the  potential  due  to  quantum  confinement. 
Due  to  the  periodicity  of  V(r),  the  solution  of  the  above  equation  are 
given  Bloch  waves  of  the  form 

'F  = eik  r u(k,r)  , (4.2) 

where  k is  the  wavevector  of  the  electron,  and  u is  a Bloch  function 
with  the  special  property  that  it  is  periodic  with  the  crystal  lattice.  In 
defining  localized  solutions  of  (4.1)  such  as  those  in  quantum  wells,  it 
is  useful  to  consider  linear  combinations  of  Bloch  wave  solutions  in 
(4.2).  Using  an  arbitrary  set  of  expansion  coefficients,  A(k),  we  can 
express  a spatially  localized  wavefunctions  as  [Zory92] 


'F  = A(k)  1 eik-ru(k,r)  d3k  = u(o,r)  J A(k)eik.rd3k  = F(r)u(r)  (4.3) 
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The  above  description  of  localized  states  is  known  as  the  envelope 
function  approximation  [Lutt55],  The  schematic  energy  band 
structure  of  a GRIN-SCH-SQW  is  shown  in  Fig.  4-1. 

For  a non-degenerate  energy  band,  such  as  the  conduction 
band,  it  has  been  shown  by  Luttinger  and  Kohn  [Lutt55]  (using  a k.p 
formalism)  that  an  effective  mass  equation  or  Schroedinger-like 
equation  for  the  envelope  function,  Fe,  as  defined  in  (4.3)  can  be 
given  as 

( Hc  + V ) = Ec  Fe,  (4.4) 


where  the  Hamiltonian  for  the  conduction  band  is  simply  given  by 


Hc  = 


-ft  2 d2 
2mc  dz2 


(4.5) 


The  Bloch  functions  have  been  removed  from  the  equation,  and  the 
effect  of  the  periodic  potential  arising  from  the  crystal  lattice  is  now 
replaced  by  a conduction  band  effective  mass,  mc.  In  this 
approximation,  the  quantum  well  created  by  the  interfacing  of  three 
materials  of  different  band  gap  truly  becomes  of  a textbook  particle 
in  a box  problem  with  F as  the  wavefunction,  and  the  material  band 
edges  as  the  potential,  V. 

The  simplicity  of  the  band  structure  in  the  conduction  band  of 
a quantum  well  relies  on  the  assumption  that  the  interaction  with 
other  energy  bands  is  weak  enough  that  we  can  treat  it 
perturbatively  by  replacing  that  interaction  with  a conduction  band 
effective  mass. 
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Fig.  4-1.  Schematic  energy  band  structure  of  a GRIN-SCH-SQW  laser. 
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However,  for  bands  degenerate  in  energy,  the  assumption  of  weak 
interaction  is  a poor  one  and  cannot  be  used.  Therefore,  the  effective 
mass  equation  (4.4)  must  be  modified  to  include  the  strong 
degenerate  band  interaction.  A modified  derivation  for  the 
degenerate  band  effective  mass  equation  has  also  been  treated  by 
Luttinger  and  Kohn  [Lutt55].  In  this  case,  we  still  obtain  an  effective 
mass  equation  for  each  degenerate  band  similiar  to  (4.4),  however,  as 
a result  of  the  degeneracy,  a coupling  term  is  introduced  which 
couples  the  equations  together.  For  the  degenerate  HH  and  LH  bands 
near  the  band  edge,  this  implies  that  we  must  work  with  four 
coupled  effective  mass  equations  (we  must  include  the  spin 
degeneracy).  Interaction  between  the  HH  and  LH  bands  is  by  far  the 
strongest,  and  we  do  not  pay  a large  penalty  by  neglecting  the  other 
bands.  The  four  coupled  effective  mass  equations  can  be  greatly 
simplifid  using  a method  suggested  by  Broido  and  Sham  [Broi85, 
Ahn88,  Ahn90].  The  two  coupled  effective  mass  equations  for  the 
degenerate  bands  can  be  expressed  as  [Zory92] 


where  Fh,i  are  the  heavy  and  light  hole  envelope  functions,  and  W+  is 
the  hermitian  conjugate  of  W.  Let  choose  the  quantum  well  direction 
to  be  along  z because  the  effective  mass  equations  would  not 
decouple  at  the  band  edge  otherwise  (see  (4.11)).  In  this  case,  kz  is 
directed  along  the  confinement  axis  and  kt  is  the  transverse  k vector 
component  and  lies  in  the  plane  of  the  well,  as  shown  in  Fig.  4-2. 


( Hh  + V)  Fh  + W Fi  = Ev  Fh 
(Hi  +V)Fi  +W+Fi=EvFi 


(4.6) 

(4.7) 
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Wave  vectors  kz  and  kt  are  two  perpendicular  k vector  components, 
with  kz  directed  along  a <100>  direction,  and  kt  directed  either  along 
a <100>  direction  or  a <110>  direction  within  the  plane  of  the  well. 
With  this,  we  can  write  [Zory92] 


The  material  constants  y 1 ,2  are  referred  to  as  the  Luttinger 
parameters  [Lutt56],  and  are  easily  related  to  the  HH  and  LH 
effective  masses,  mhh  and  mih  through  (4.10).  A third  Luttinger 
parameter,  73,  exists  in  the  coupling  term,  W.  The  coupling  term 
takes  a slightly  different  form  when  kt  is  directed  along  {100}  and 
{110}  directions  [Zory92}.  We  can  define  the  two  forms  as 


Because  the  energy  bands  are  different  in  the  {110}  and  {100}  plane, 
and  because  electrons  will  of  course  exist  with  k vectors  in  all  planes, 
we  use  here  the  axial  approximation  [Alta85],  which  assumes  72  = 73 
only  in  the  W term. 


Hh  = (7l  “ 272)  kz2  + (7!  + 272)  kt2, 
Hi  = (71  + 272)  kz2  + (71  - 272)  kt2, 


(4.8) 

(4.9) 


where 


(4.10) 


W = V3  kt  (72  kt  - i2  73  kz)  for  {100}  planes,  (4.11a) 

W = V3  kt  (73  kt  - i2  73  kz)  for  {110}  planes.  (4.11b) 
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<100> 


<100> 


Fig.  4-2.  Coordinate  system  used  in  the  valence  band  model.  The 
transverse  or  in-plane  k vector,  kt,  can  be  directed  along 
either  a <100>  or  <11 0>  direction,  whereas  the  confinement 
axis  must  be  along  a <100>  direction. 
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From  Eqs.  (4.8)-(4.1 1),  we  have  two  coupled  differential 
equations  for  Fh(z)  and  Fi(z)  given  by 

[ (.2  (Yl  + Y2)  kt2  - 2 (Yl  - 272)  + V(z)  ] Fh(z) 

+ (y  7 kt2  - t/3  73  kt  F|(z)  =^E,  Fh(z),  (4.12a) 

(-J-  7 ki2  + 'll  73  kt ")  Fh(z)  + [ (7i  - 72)  kt2 

- J <71  + 272)  ^ ) + ^ VW  ) Fl(z)  = ^ Ev  Fi(z),  (4.12b) 


It  is  convenient  to  normalize  Eqs.(4.12a)  and  (4.12b)  by  defining 


H}.ql2 

H2 


Ev, 


r_  m0L2 
f-  U2 


V($), 


Ai  = 2 (Yl  + 72)  L2kt2,  A2  = ^ (Yl  - Y2), 
C=-y-  Y L2  kt2,  D=  V3  Y3  Lkt , Y = ^ (72  + 73) 
Bi  = 2 (Yl  - 72)  L2kt2,  and  B2  = J (Yl  + Y2)- 


(4.13) 


Then  Eqs. (13. a)  and  (13. b)  become 

(Ai  - A2  ^ + f(«)  Fh(0  + (C  - D ~ ) F,(4) 

= hi  Fh(^)  (4.14a) 

(C  + D ^ ) Fhft)  + (Bi  - B2  ^ + ffc)  ) Fi($) 


= Xi  Fj(^) 


(4.14b) 
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By  the  same  procedures,  conduction  band  effective  mass  equation 
(4.4)  becomes 


(- 


J_  ^2_ 

mc  d^2 


+ 


2m0L2 

H2 


■v«)  + 


kt2L2 

mc 


) Fc(0  = Ec  Fc(4) 


(4.15) 


where  mc=  mc*/m0,  and  Ec  is  the  conduction  subband  energy. 

4.2.2.  Transition  matrix  element. 

The  transition  of  the  optical  matrix  elements  depend  on  the 
polarization  of  the  incident  light  and  we  consider  two  cases:  TE  and 
TM.  The  transition  matrix  element  can  be  simplified  in  terms  of 
normalized  envelope  wavefunction 

4 

IMjl2  = lM0l2  [ I <Fi  I Fc>l2  ] TMpolarization  (4.16) 

IMTI2  = IM0I2  [ l<Fh  I Fc>l2  + j l<Fi  I Fc>l2  ] TE  polarization  (4.17) 

where  IM  0 12  defined  in  Eq.(3-8)  is  the  bulk  material  transition  matrix 
element.  Here  we  used  m*  = 0.053mo  in  lMQl2for  QW  [Yan90].  In 
making  the  envelope  functions  and  evaluating  the  overlap  integrals 
numerically,  we  must  make  sure  that  they  are  properly  normalized. 
Normalization  of  the  wavefunctions  is  obtained  through  the  following 
relations 

p. 

Fj  (norm)  = where  Nh,i  = < Fhl  Fh>  + < Fj  I Fj  >, 


Ne  = < Fe  I Fe  > . 


(4.18) 
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4.2.3.  Expression  for  Gain  and  Spontaneous  Emission 
4. 2. 3.1.  Gain 

The  correct  expression  for  the  gain  for  bulk  or  QW  material  can 
be  written  in  MKS  units  as  [Yan90] 

g(E)  =[Cjj-  ) (—  )2  (-  4co2  ) I MT  I 2 Prd(E)  ] (4.19) 

.[f  ] [fc(l-  fv)-fv(l-  fc)] 

=(£^i)lMT'2^(E)(fc-w  (4-2o> 

The  first  term  in  (4.19)  is  an  expression  for  Fermi's  Golden  Rule 
which  expresses  the  probability  per  unit  time  for  transitions  occuring 
from  a filled  state  to  an  empty  state.  The  second  term  converts  the 
probability  per  unit  time  into  probability  per  unit  length.  The  third 
term  represents  the  probability  that  the  initial  state  is  filled  and  the 
final  state  is  empty.  IMjl2  is  the  transition  matrix  element  in  QW  and 
Prd(E)  is  the  reduced  density  of  states,  which  is  given  by  [Made78] 

p(kt) 

prd(E)=  [dEeh/  dkt  ] Eeh=hco  (4.21) 

where  p(kt)  is  density  of  states  and  Eeh  is  the  transition  energy.  For 
the  parabolic  energy  band,  we  can  get  the  closed  form  of  reduced 
density  of  states.  For  the  nonparabolic  bands  as  in  the  valence  band 
of  quantum  wells,  prcj  cannot  be  given  in  closed-form.  The  electron 
and  hole  occupation  probability,  which  is  determined  from  the 
Fermi-Dirac  distribution  for  electrons  in  the  conduction  and  valence 
bands,  is  defined  as 
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fc’v  - 1 + exp  [ ( E e . h - Efc,fv)/kBT]  ' 


(4.22) 


where  EfC)fv  are  the  non-equilibrium  quasi-Fermi  levels  in  the 
conduction  and  valence  bands.  For  the  parabolic  energy  band, 
transition  energy  can  be  written  as 


, -&2k2 

Eeh  = Ee  - Eh  = Eg  + 


ti2k2 

+ 2mv  ’ 


(4.23) 


and  the  individual  energies  Ee,h  are  related  to  the  transition  energy, 
Eeh,  by 


Ee  = Ec  + (Eeh  - Eg')  ™r  , Eh  = Ev  - (Eeh  - Eg'V^- 

ill Q III  Y 


(4.24) 


The  primed  band  gap  energy,  Eg'  = Ec  - Ev,  is  defined  as  the  band  gap 
between  two  given  subbands  in  a QW  and  is  just  the  material  band 
edge  in  bulk  material.  The  reduced  mass,  mr,  is  given  by 


1 _ 1 1 
m r “ mc  + m v 


(4.25) 


4. 2. 3. 2.  Spontaneous  Emission 

The  expression  for  the  spontaneous  emission  rate  is  useful  in 
calculating  the  current  required  to  provide  a given  gain  found  in 
(4.20).  For  the  case  of  k-selection,  one  must  find  the  radiative 
component  of  the  current  density  by  evaluating  the  spontaneous 
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emission  rate  integrated  over  all  energies.  The  spontaneous  emission 
rate  per  unit  energy  at  a given  energy  ( rsp(E)  has  units  of  per  unit 
time  per  unit  volume  per  unit  energy)  is  given  by  [Yan90] 

rsp(E)  = (f)(-)2  (-^i“  ) I Wave  I 2 Prd(E) 

,D(E).  fc(l  - f„)  (4.26) 

where  D(E)  is  the  optical  mode  density,  in  the  material  with  a 
refractive  index  of  n,  given  by 

n3E2 

D(E)=  7t2  1i 3 c3  (4'27) 

, and  iMavel2  is  the  average  momentum  matrix  element,  defined  by 
[Zory92] 


I Mave  1 2 = J X |MtI2.  (4.28) 

all  three 
polarization 

In  QW  structures,  the  matrix  element  is  enhanced  for  certain  electric 
field  polarizations  and  reduced  for  others  [Yama84].  For  total 
spontaneous  emission  we  are  interested  in  the  total  output  of  light, 
not  which  polarizations  of  light  are  being  emitted  spontaneously. 
Thus,  we  take  the  average  strength  of  the  transition  matrix  element 
over  all  three  polarization.  The  integrated  spontaneous  emission 
output  is  just  equal  to  the  total  radiative  recombination  of  carriers 
and  is  therefore  equal  to  the  radiative  component  of  the  current 
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density  required  to  achieve  a quasi-Fermi  level  separation.  The 
relation  can  be  written  as 


where  d is  the  QW  thickness. 

4. 2. 3. 3.  Carrier  density  and  Quasi-Fermi  level 

The  carrier  density  in  a given  band  can  be  found  for  a given 
quasi-Fermi  level  by  integrating  the  density  of  states  multiplied  by 
the  occupation  probability  over  the  entire  band.  For  the  parabolic 
subbands,  such  as  in  the  conduction  band  of  a QW,  we  can  write  the 
electron  concentration  as 


The  sum  is  over  all  quantized  subbands  within  the  conduction  band 
of  the  quantum  well,  and  the  Ec  are  the  quantized  energy  levels. 

For  non-parabolic  bands,  such  as  in  the  valence  band  of  a QW, 
where  subband  structures  are  nonparabolic,  Eq.(4.30)  is  no  longer 
valid  because  density  of  states  becomes  a complicated  function.  In 
this  case,  it  is  more  appropriate  to  find  the  carrier  density  by 
numerically  integrating  over  k-space,  because  k states  are  always 
uniformly  distributed  in  k-space,  independent  of  how  complicated 
the  energy  bands  become.  For  this  case,  we  have  [Zory92] 


(4.29) 


N = ^^fTd  £ In  { 1 + exp[  -(Ec  - Efc)/kBT  ]} 


(4.30) 


n 
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kmax 

P = 2 S Jp(kt)  [1  - fv(Ev(k,))]  dk,  (4.31) 

0 

n 

where  p(kt)  = is  density  of  states.  The  sum  is  again  over  all 

quantized  subbands,  Ev(kt),  within  the  valence  band  of  the  QW, 
including  both  heavy  and  light  hole  bands.  kt  is  the  k vector 
component  in  the  plane  of  the  well  and  kmax  represents  some 
numerical  limit  beyond  which  the  contribution  to  the  integral  can  be 
neglected. 


4.3.  Numerical  procedures  and  results 
The  multiband  effective  mass  equations  were  solved  numerically 
using  the  finite-difference  method.  This  method  is  very  convenient 
in  that  an  arbitrary  QW  geometry  or  potential  profile  can  be  handled 
easily,  but  the  drawback  is  that  the  size  of  the  matrix  becomes  larger 
for  numerical  accuracy.  After  discretizing  the  differential  equations 
and  evaluating  the  coefficients  in  (4.14a),  (4.14b),  and  (4.15)  for  each 
kt,  the  problem  reduces  to  an  eigenvalue  problem  to  be  solved  for 
the  subband  energies  and  the  the  envelope  wavefunction.  The 
material  parameters  to  be  used  in  the  calculations  are  given  in  Table 
III.  Most  bulk  parameters  were  taken  from  reference  [Land82],  with 
the  exception  of  the  AlGaAs  bandgap  which  was  taken  from  Casey 
and  Panish  [Case78].  The  conduction  band  offset  in  the  GaAs/AlGaAs 
system  was  assumed  to  be  62%,  taken  from  Watanabe  [Wata85].  For 
the  Luttinger  parameters  of  GaAs  we  used  the  value  of  yi  = 6.85,  72  = 
2.1,  73  = 2.9.  The  effective  masses  of  electrons  and  holes  in  the 
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barrier  material  were  assumed  to  be  the  same  as  those  in  the  QW 
because  of  the  conversion  problem.  The  difference  on  the  Luttinger 
parameters  and  electron  effective  masses  in  QW  and  GRIN  regions 
may  introduce  small  corrections  to  the  wavefunctions  and  to  the 
subband  energy  levels  of  the  order  of  a few  meV,  which  can  be 
neglected  as  compared  with  the  band-gap  energy  [Sand87].  Boundary 
value  ( z = ± 46.875  nm  ) should  be  chosen  such  that  node  point 
could  be  matched  at  GRIN/QW  interface  to  make  the  potential  barrier 
change  smooth  and  eigenfunction  go  to  zero  at  boundary.  At 
GRIN/QW  interface  we  used  averaged  potential  value  of  the 
GRIN/QW  interface.  Step  size(=dz)  was  also  critical  for  the  numerical 
conversion  and  accuracy.  In  our  calculation,  matrix  size  was  400  x 

400  and  dz  was  0.46875  nm. 

The  E-k  dispersion  curves  for  the  conduction  and  valence 
subbands  of  7.5  nm  GRIN-SCH-SQW  calculated  by  the  k.p  techniques 
are  given  in  Fig.  4-3  and  4-4  respectively.  The  valence  subbands 

calculated  by  the  k.p  method  are  totally  nonparabolic  due  to  the 

effect  of  the  band  mixing  which  is  taken  into  account  in  this  method. 
At  kt=0,  heavy  and  light  holes  are  exactly  decoupled.  For  a nonzero  kt 
the  increasing  admixture  of  heavy-  and  light  hole  states  gives  rise  to 
strong  nonparabolicities  in  the  valence  band  structure.  Some  of  the 
bands  are  seen  to  have  negative  zone  center  effective  masses. 
Subbands  are  labelled  as  heavy-hole  or  light-hole  (Fig.  4-5) 

according  to  their  envelope  wave  function  shape  at  band  edge  where 
in-plane  wavector  kt  is  equal  to  zero. 
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Table.  Ill 

AlxGai-xAs  bulk  parameters 

Eg  (eV)  = 1.424  + 1.247x 

Eg  (eV)  = 1.424  + 1.247x  +1.147(x-0.45)2 

mc  = 0.067  + 0.083x 


x < 0.45 
x > 0.45 


Yl  = 6.85  (1  -x  ) + 3.45  x 

72  = 2.10  (1-x  ) + 0.68  x 

73  = 2.90  (1-x  ) + 1.29  x 


(meV) 
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k (1/A) 


Fig.  4-3.  Conduction  subband  structure  of  a 7.5  nm  GaAs/AlGaAs 
GRIN-SCH-SQW. 
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Fig.  4-4.  Valence  subband  structure  of  a 7.5  nm  Ga/AlGaAs  GRIN- 
SCH-SQW. 


Amplitude 
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Fig.  4-5.  Hole  subband  envelope  wave  function  at  band  edge. 
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In  section  4.2.2,  we  discussed  polarization  dependent  effects  in 
quantum  wells.  We  provided  a general  method  for  predicting  band 
edge  transition  strengths.  We  can  find  the  TM  and  TE  transition 

strengths  from  Eqs.(4.16)  and  (4.17)  for  each  subband  transition  as  a 
function  of  in  plane  k vector.  Figs.  4-6,  4-7,  4-8,  and  4-9  show  TE  and 
TM  transition  strength  associated  with  the  energy  eigenfunctions 
between  conduction  and  valence  subbands.  These  figures  show  the 
selection  rules  for  different  polarizations  of  the  optical  radiation. 
Only  at  kt  =0  is  the  An  = 0 selection  rule  seen  to  hold  and  it  is 

important  to  note  that  all  of  the  band-to-band  transitions  have 
nonzero  matrix  elements.  The  forbidden  transitions  are  allowed  in 
both  mode  transition  due  to  the  redistribution  of  the  transition 

strength  among  the  various  subband  transitions.  The  optical  matrix 
elements  are  seem  to  be  strongly  kt  dependent  due  to  valence  band 
mixing. 

The  Fermi  levels  found  for  a given  carrier  density,  N or  P,  by 
assuming  charge  neutrality,  N=P,  are  shown  in  Figs.  4-10  and  4-11. 
Using  the  Eqs.(4.30)  and  (4.31),  the  conduction  and  valence  band 

quasi-Fermi  energy  level  Efc  and  Efv  with  respect  to  electron  (N)  and 
hole  carrier  density  (P)  was  calculated  by  integrating  the  charges 
within  the  conduction  and  valence  bands  numerically.  The  density  of 
states  can  be  calculated  from  the  E-k  diagrams  of  Fig.  4-3  and  4-4. 
The  values  of  the  density  of  states  without  including  spin  degeneracy 
are  defined  by  p(kt)  = . To  simplify  the  final  gain  calculation,  we 

used  curvefitted  equations  for  the  resulting  Efc  vs  N and  Efv  vs  P 
using  polynomial  degree  of  8.  IMSL  library  PPITT  was  used  for  the 


Transition  Strength 


71 


k (1/A) 


Fig.  4-6.  Transition  strength  for  TE  light  polarization  between  1st 
conduction  and  hole  subbands. 
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Fig.  4-7.  Transition  strength  for  TE  light  polarization  between  2nd 
conduction  and  hole  subbands. 
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Fig.  4-8.  Transition  strength  for  TM  light  polarization  between  1st 
conduction  and  hole  subbands. 
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Fig.  4-9.  Transition  strength  for  TM  light  polarization  between  2nd 
* conduction  and  hole  subbands. 
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Fig.  4-10.  Conduction  band  quasi-Fermi  energy  level,  Efc  vs 
electron  carrier  density,  N 
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Fig.  4-11.  Valence  band  quasi-Fermi  energy  level,  Efv  vs  hole  carrier 
density,  P 
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numerical  integration  of  carrier  density.  However,  we  could  not  use 
the  above  curve  fitted  equations  for  the  estimation  of  reduced 
density  of  states  because  derivative  of  E with  respect  to  k produced 
big  error.  We  used  cubic  spline  method  (IMSL  library  PPDER)  which 
gave  exact  value  of  reduced  density  of  states. 

We  then  calcaulated  the  theoretical  gain  using  standard 
methods  [Yan90,  Zory92].  Since  some  of  the  subbands  are  double 
valued  in  energy,  the  gain  computations  are  done  in  the  k space. 
Therefore,  the  gain  is  first  found  between  conduction  and  valence 
subbands  in  the  k space  within  the  QW  plane,  and  then  summed  for 
the  transitions  having  the  same  energy  by  using 


where  g'(kt)  is  the  gain  for  each  kt  increment  which  is  given  in 
Eq.(4.20),  and  d(E'-E)  is  the  delta  function.  The  g(E)  is  then  the 
summed  gain  spectra  for  each  pair  of  electron  and  hole  subbands. 
The  total  gain  is  then  found  by  simply  adding  such  gain  spectra  for 
each  pair  of  subbands  (IMSL  library  PPDER).  The  results  of  the  TE 
peak  gain  as  a function  of  energy  are  shown  in  Fig.  4-12  and  4-13  for 
two  different  carrier  density.  These  figures  show  that  the  forbidden 
transitions  are  partially  allowed  due  to  band  mixing,  and  the 
transition  shift  takes  place  to  n=2  subbands  at  high  carrier  density. 


OO 


(4.32) 
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Fig.  4-12.  TE  mode  gain  as  a function  of  energy  at  carrier  density 
3xl018  cm-3. 
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Fig.  4-13.  TE  mode  gain  as  a function  of  energy  at  carrier  density 
3xl019  cm*3.  . 
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The  results  of  the  peak  gain  vs  J/d  calcualation  at  25°C  for  the 
TE  mode  are  shown  in  Fig.  4-14  assuming  all  the  injected  carriers 
recombine  radiatively  (ijsp=l)  (black  circle:  experimental  data).  For 
comparison  purposes,  we  have  included  parabolic  band 
approximation  by  Chinn.  Parabolic  model  predicts  a radiative  current 
density  similar  to  the  valence  band  mixing  model  for  low  gain 
regions.  However,  valence  band  mixing  model  needs  more  current  in 
high  gain  region  than  by  the  parabolic  band  model.  Fig.  4-15  shows  g 
vs  J/d  for  the  parabolic  and  nonparabolic  band  strucutre  by  taking 
riSp  into  account.  We  used  measured  injection  efficiency  r|jn=0.82  for 
this  calculation  [Zmud91].  For  Auger  recombination  rate  calculation, 
we  used  same  formula  as  in  Chinn  model  [Chin88].  In  Fig.  4-16  we 
compared  g vs  J/d  at  125°C  with  that  of  25°C  assuming  rjsp=l. 
Experimental  data  is  for  L=500  pm  laser  and  we  used  oti  value 
measured  at  room  temperature.  We  were  unable  to  estimate  oci  and 
rj  in  at  this  high  temperature  because  Pd  vs  L measurements  were  not 
linear. 


4.5.  Summary 

Combination  of  the  effects  by  valence  band  mixing  and 
measured  injection  efficiency  fit  the  7.5  nm  QW  experimental  data, 
leaving  no  room  for  the  CS  contribution.  Recent  theoretical  studies 
shows  that  the  CS  contribution  to  the  threshold  current  is  not  as  big 
as  reported  earlier  [Yama81].  It  is  found  that  the  line  shape  functions 
in  low  dimensional  systems  have  a very  small  homogeneous 
broadened  widths  and  strong  convergent  characteristic  [Yama87, 
Asad89,  Ohto91].  Leakage  current  density  estimated  in  terms  of 
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measured  injection  efficiency  (riin)  is  significant,  comparable  to  the 
current  density  contribution  by  carrier  scattering  in  broadening 
model. 
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Fig.  4-14.  Peak  TE  mode  gain  vs  J/d  assuming  all  the  injected 
carriers  recombine  radiatively  (Tisp=l)  (black  circle: 
experimental  data) 
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Fig.  4-15.  Peak  TE  gain  vs  J/d  for  the  parabolic  and  nonparabolic 
band  strucutre  by  taking  rj  sp  into  account. 
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Fig.  4-16.  Threshold  gain  as  a function  of  J/d  with  r|sp  = 1 at  125°C. 


CHAPTER  5 

SUMMARY  AND  CONCLUSIONS 


The  goal  of  this  study  was  to  integrate  the  results  of 
experimental  and  theoretical  work  to  come  up  with  a model  capable 
of  predicting  threshold  current  density  (Jth)  at  any  temperature  for 
QW  lasers.  As  discussed  in  Chapter  1,  it  was  necessary  to  do  this 
since  the  existing  CZR  model  for  QW  lasers  [Chin88,  Chin89]  was 
found  to  be  deficient. 

To  sort  out  the  problem,  we  chose  two  cases,  thick  40nm  well 
(bulk)  and  thin  7.5  nm  QW.  In  thick  wells,  all  the  complicated 
phenomena  related  to  the  QW  vanish  and  only  CS  is  expected  to  be 
fundamental.  We  investigated  the  effect  of  CS  on  Jth  using  standard 
bulk  theory.  In  thin  7.5  nm  QW  where  quantum  size  and 
nonradiative  effects  are  important,  we  studied  the  effect  of  CS  on  Jth 
using  new  VBM  QW  theory.  The  leakage  current  was  included  into 
the  total  current  in  terms  of  spontaneous  emission  efficiencies  (r)sp), 
which  can  be  expressed  as  the  product  of  two  efficiencies,  r|in  = 
injection  efficiency  and  qr  = active  layer  radiative  efficiency. 
Injection  efficiency  was  estimated  from  the  differential  quantum 
efficiency  versus  laser  cavity  length  measurements.  Radiative 
recombination  rate  was  calculated  theoretically. 

From  the  above  work  we  found  that  a reasonable  fit  to  the 
thick  QW(bulk)  experimental  data  at  both  25°C  and  125°C  could  be 
obtained  using  strict  k-selection  bulk  theory  and  a knowledge  of  the 
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carrier  injection  efficiency(r|in).  In  making  the  fit,  carrier  scattering 
(CS)  effect  is  not  necessary,  quite  contrary  to  the  Yamada's  results 
[Yama81].  The  essence  of  the  problem  appears  to  be  that  if  measured 
J th  values  are  low  enough,  Stern’s  theory  without  CS  is  sufficient  to 
explain  results.  This  finding  suggested  the  possibility  that  the  effect 
of  CS  on  Jth  might  not  be  important  even  in  thin  QW  lasers.  This  was 
verified  with  VBM  QW  theory  combined  with  a knowledge  of  the 
carrier  injection  efficiency(ri  jn).  Combination  of  the  effects  by 
valence  band  mixing  and  measured  injection  efficiency  fit  the  7.5  nm 
QW  experimental  data,  leaving  no  room  for  the  CS  contribution. 
Leakage  current  density  estimated  in  terms  of  measured  injection 
efficiency  (rjin)  is  significant,  comparable  to  the  current  density 
contribution  by  CS  in  the  CZR  model.  This  does  not  necessarily  rule 
out  the  possibility  that  experimental  spectra  can,  in  principle,  be 
obtained  which  are  subject  to  the  strict  k-selection  rule.  However, 
recent  theoretical  studies  shows  that  the  CS  contribution  to  the 
threshold  current  is  not  as  big  as  reported  earlier  [Yama81].  It  is 
found  that  the  line  shape  functions  in  low  dimensional  systems  have 
a very  small  homogeneous  broadened  widths  and  strong  convergent 
characteristic  [Yama87,  Asad89,  Ohto91]. 

It  is  concluded  that  good  predictions  of  threshold  current 
density  at  any  temperature  can  be  made  for  QW  lasers  provided  one 
uses  QW  theory  developed  in  this  work  and  one  knows  how  to 
predict  the  temperature  dependence  of  carrier  injection 
efficiency (ri in).  In  our  low  temperature  QW  experiments,  we  were 
able  to  determine  Tiin  from  the  differential  quantum  efficiency  (rj d)  vs 
L measurements  and  get  good  agreement  between  theory  and 
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temperature.  However,  we  were  unable  to  estimate  r\  i n at  high 
temperature  because  rid  vs  L measurements  were  not  linear 
indicating  that  r\  in  was  a function  of  cavity  length  at  high 
temperature.  As  a consequence,  r|in  at  high  temperature  for  the 
30/70  QW  structures  studied  here  remains  to  be  determined.  We  are 
of  the  opinion  that  the  correct  technique  for  calculating  rj  ± n has  been 
published  recently  by  Hirayama  et  al.  [Hira92].  It  will  be  interesting 
to  see  if  the  combined  theories  (the  VBM  theory  developed  in  this 
work  and  that  of  Harayama  et  al  for  pjn)  are  able  to  make  correct 
predictions  of  Jth  for  QW  structures  at  any  temperature. 
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